The optical phonon spectrum of SmFeAsO 
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We measured the Raman and the Infrared phonon spectrum of SmFeAsO poly crystalline sam- 
ples. We also performed Density Functional Theory calculations within the pseudopotential ap- 
proximation to obtain the structural and dynamical lattice properties of both the SmFeAsO and the 
prototype LaFeAsO compounds. The measured Raman and Infrared phonon frequencies are well 
predicted by the optical phonon frequencies computed at the V point, showing the capability of the 
employed ab-initio methods to describe the dynamical properties of these materials. A comparison 
among the phonon frequencies of different oxypnictides suggests a possible role of the high frequency 
phonons in the pairing mechanism leading to superconductivity in these materials. 
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I. INTRODUCTION 

The recent discovery of high temperature superconductivity in iron-based oxypnictides [l[ has triggered in the 
scientific community a strong research effort on this new class of superconductors. Superconducting pnictides are 
layered compounds containing a transition metal (Fe, Co, Ni) bonded with a group V element (As and P) and 
intercalated with layers formed by a rare-earth (RE = La, Sm, Nd, Ce) and oxygen. Superconductivity appears upon 
doping with electrons (substitution of O with F) or with holes (oxygen deficiency), while the parent compounds show 
an antiferromagnetic metallic ground state. After the discovery of superconductivity at 26 K in F-doped LaFeAsO 
higher critical temperatures were soon obtained bysubstituting La with Sm 0, Ce 0], Pr Q and Nd Also 
oxygen- free compounds like ^4Fe2As2 (A = Ca, Sr, Ba [aj7|) were found to exhibit superconductivity, even if at lower 
temperatures, thus opening the way to a new 'iron-age' [8|. 

The prototype LaFeAsO shows an anomaly near Tsdw = 150 K both in the resistivity and in the magnetic 
susceptibility, which is attributed to a spin-density- wave (SDW) instability [H [HI El . Just above Tsdw, the system 
under goes a structural transition from a high temperature tetragonal (P4/nmm) phase to an orthorhombic (Cmma) 
phase |12l.[l3j. It is generally re cog nized that doped compounds become superconducting when electron or hole doping 
suppresses the SDW instability [lj]. However, a satisfactory theory of the superconducting mechanism in the pnictides 
has not yet been found. In particular, the role of phonons in the pairing mechanism has been strongly questioned due 
to the small value of the electron-phonon coupling[15]. The scenario is also complicated by difficulties in an ab-initio 
description of the normal state properties of these compounds [l6j . A careful comparison between experimental and 
calculated properties is thus timely and important in the validation of the present first-principle theoretical methods. 
Structural and dynamical properties of the i?i?FeAsO compounds are an obvious ground where this comparison is 
due. 

In this context we measured both the Raman and the Infrared (IR) phonon spectrum of pure SmFeAsO poly- 
crystalline samples, well characterized in terms of crystal structure (P4/nmm at high temperature as in SmFeAsO), 
magnetic, thermal and transport properties showing that Tsdw ~ 135-^140 K [l?], EH- We have then calculated the 
phonon spectrum and compared the theoretical results with those of the experiments. 
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II. RAMAN AND INFRARED SPECTRA 



Raman spectra were measured in back-scattering geometry using a confocal micro- Raman spectrometer equipped 
with a He-Ne laser at 632.8 nm. Measurements were performed on fresh-cut surfaces without any additional treatment. 
A camera coupled to the microscope allowed to select the point over the sample surface on which to focus the laser 
spot (about 5 /im by using a 20X objective). Since SmFeAsO is a very poor Raman scatterer fl9| . spectra with a 
good signal-to-noise ratio were collected within the 80-1000 cm' 1 range by setting long acquisition times and a large 
number of spectra to be averaged. It was verified that the obtained results are not affected by laser-induced sample 
heating. The polycrystalline nature of the sample prevent any polarization analysis of the collected spectra. 

A typical room temperature Raman spectrum is shown in the inset of Fig[T] where Raman peaks are evident between 
170 and 210 cm -1 . Only in a few cases the collected spectra show Raman contributions due to samarium oxide [201 ] . 
in agreement with the results of the analyses performed on the same sample [13] • Spectra collected from different 
points (P1-P4) are reported in Fig. [T]in a limited frequency range to highlight the details of the Raman peaks. Three 
lines are evident, at frequencies (170, 201 and 208 cm -1 ) which coincide with those reported in a previous room 
temperature Raman study which includes polarization analysis [19]. The different relative intensities of the peaks we 
observed in different points must be ascribed to polarization effects. Indeed the intensity of the Raman peaks can 
depend on the angle between the laser polarization and the possible preferred orientation of the crystallites in the 
scattering volume. 

We also studied the temperature dependence of the Raman spectrum between 300 and 100 K by mounting the 
sample in a cryostat, which requires longer acquisition times. The comparison between the spectra collected from the 
same point (P4) at 300 and 100 K (see Fig. [I} shows no remarkable changes in the number of phonons or in the peak 
frequencies occurs decreasing temperature below Tsdw, besides the expected phonon hardening and narrowing. The 
300 and 100 K peak frequency values are reported in Table [V] 

The infrared reflectivity R(w) of the same SmFeAsO sample was then measured from 50 to 5000 cm -1 and between 
17 and 300 K. The reference was the sample itself, gold-coated by evaporation in-situ. These measurements, performed 
using a Michelson interferometer at the SISSI Infrared beamline of the ELETTRA Storage Ring (Trieste, Italy), showed 
that R(w) is independent of temperature above 1000 cm . The reflectivity measurements were then extended at up 
to 30000 cm" 1 at room temperature, by using a monochromator and an aluminum film as reference. 

The inset of Fig. [2] shows that, as expected, R(w) has a bad metallic character. A Hagen- Rubens relation was then 
employed for the extrapolation to u = 0, usual power laws for that to lu — oo, and the optical conductivity ai(u) 
was then obtained by standard Kramers-Kronig transformations. Due to the wide frequency range where R(w) is 
measured directly, (J\(uj) in the phonon frequency range (see Fig. [5]) does not depend appreciably on the choice of the 
extrapolations. 

Although the present paper is devoted to the study of the optical phonons, it is worth noting that, in agreement 
with previous results on LaFeAsO @, the optical conductivity of SmFeAsO does not provide a clear evidence of a 
SDW gap, indicating that not all of the Fermi surface is gapped. Only an increase of the low frequency conductivity 
owing to the increase of carrier mobility is observed for decreasing temperature below Tsdw [HI- 

In the far-infrared spectrum five phonon peaks at 101, 259, 276, 377, and 448 cm" 1 are evident at low temperatures. 
Above 150 K, the 259 and 375 cm" 1 peaks are hardly discernible and a standard Lorentzian fit was employed to 
determine their peak frequencies. As in the Raman case, no remarkable change in the peak frequencies occurs on 
decreasing temperature. The 300 and 100 K peak frequencies are reported in Table [V] We note that the structural 
phase transition around Tsdw from tetragonal to orthorhombic recently observed in SmFeAsO [22], as well as in 
LaFeAsO [H Il3| and in NdFeAsO[23|], does not significantly affect the phonon frequencies, as also observed in the 
LaFeAsO case[9j,ll0(. However, the significant increase in the relative intensities of the 259 and 377 cm -1 peaks below 
150 K may be regarded as an effect of the structural transition around Tsdw- 



III. CALCULATIONS 



The experimental data reported in Table [V] were compared with Density Functional Theory (DFT) calcula- 
tions of the SmFeAsO optical phonons, performed within the pseudopotential approximation for the electron-ion 
interaction[U[2^, [2(|. The f-states of Sm were considered to be inside the core and a formal valence +3 was assumed 
for Sm (6s 2 5c? 1 as valence electrons, with 5s and 5p treated as semi-core states). The electronic density, integrated on 
a uniform mesh of 14 x 14 x 10 k-points, converged within an energy cutoff of 400 eV. A local-density approximation 
with gradient corrections 27j was used for the exchange correlation energy. 

For comparison and cross checking with previous linear-response theory calculations, we also performed calculations 
for the prototype LaFeAsO compound, in the same P4/nmm structure as SmFeAsO. The theoretical lattice constants 
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were obtained by minimizing the total energy as a function of the volume, after relaxing all the internal degrees of 
freedom. In all calculations, both systems are considered in the paramagnetic phase. 

The resulting lattice constants for LaFeAsO are a = 4.017 A, c = 8.677 A, z^ a = 0.144 and za s = 0.638, in good 
agreement with previous computations [l5j], while for SmFcAsO we find a = 3.948 A, c = 8.336 A, z Sm = 0.141 and 
zas = 0.647, in reasonable agreement with the experimental findings [17| (for a discussion on the agreement between 
experimental and computed parameters see Ref [lllHsl])- We notice that the smaller Sm ion reduces both the lattice 
constants: the in-plane one is reduced by about 2%, while the out-of-plane one is reduced by 4%. This reduction 
does not reflect on all interatomic bond lengths. The most affected parameter is the height of RE atoms above the 
O plane, passing from 1.25 A in LaFeAsO to 1.17 A in SmFeAsO. Consequently, the RE-0 bond length is reduced 
from 2.36 A to 2.30 A. The Fe-As bond lenghts instead changes from 2.34 A in LaFeAsO to 2.36 A in SmFeAsO, 
although the As height above the Fe plane increases upon Sm substitution. As we will see, these structural properties 
are crucial to explain the differences in the phonon frequencies between these compounds. 

The phonon frequencies at T, to be compared with the measured frequencies, were obtained by diagonalizing 
the dynamical matrix calculated in the frozen phonon approach [30l|. For each phonon mode we report in Table IVl 
experimental and computed frequency, the irreducible representation associated with the mode, the polarization (X,Y: 
in-plane, Z: out-of-plane), the optical activity (IR or Raman), and the dominant atoms involved. A comparison with 
the present calculations allows one to assign all the observedphonon lines. We also report the calculated frequencies 
in LaFeAsO, in good agreement with previous computations [Iol.[l5l]. and the corresponding experimental values from 
literature data [3, [l!| . 

A close inspection of data in Table [V] shows a good agreement (very good in some cases) between the measured and 
calculated frequencies. As expected, the low-frequency part of the spectrum (up to about 210 cm -1 ) is dominated 
by vibrations mainly involving the massive atoms (RE and As). The comparison with the calculated frequencies 
in LaFeAsO reveals a small but systematic increase of the frequencies associated with in-plane motions. The larger 
Sm mass with respect to La would account for a slight reduction of the frequencies (u>sm/^La — (-^La/Msm) 5 ) in 
passing from La to Sm, but this is washed out by the reduction of the strong 2D in-plane bonds. The same happens 
for the mode with Z-polarization (97 cm -1 ) that also involves the Fe-As motion. On the other hand, the modes at 
177 and 203 cm -1 (out-of-plane vibrations of RE and As) turn out to be softer than the respective ones in LaFeAsO, 
as partially confirmed by the experiments. 

It is worth noting that, due to the weakness of the RE- As bond, the ionic interactions between the layers are weaker 
than the in-layer ones. We then expect that a simple rescaling of the frequencies with the RE mass should be valid. 
This was verified by using the force constants calculated in the case of LaFeAsO (at its own lattice constants) and 
diagonalizing the dynamical matrix with the La atomic mass substituting that of Sm. We found that the softening 
observed in SmFeAsO is exactly predicted by this model calculation for only these two modes. 

At frequencies up to 300 cm -1 the observed modes only involve the Fe-As network. The mode at 224 cm -1 involves 
only the Fe atoms moving along the Z-direction, and thus is expected to be unaffected by the Sm substitution and by 
the slight Fe-As bond reduction. The other modes, Fe-As stretching modes, are consistently higher in passing from 
LaFeAsO to SmFeAsO. 

Only in the region above 300 cm , which is dominated by the vibrations involving the oxygen atoms, the phonon 
frequencies of the two compounds are significantly different. In fact, the substitution of La with Sm produces a 
shortening of the RE-0 bond length and of the zre parameter (see the discussion above). This produces an overall 
stiffening of the RE-0 modes; in particular, the highest frequency E g phonon is at 434 cm -1 in LaFeAsO, at 503 
cm -1 in SmFeAsO. 



IV. CONCLUSION 

We performed Raman and IR measurement on SmFeAsO, thus obtaining accurate information on the optical phonon 
frequencies at the T point. They show that the structural phase transition around 140 K does not significantly affect 
the phonon frequencies, as previously observed in the LaFeAsO case. We also performed a theoretical study of the 
structural and dynamical lattice properties by means of accurate Density Functional Theory calculations within the 
pseudopotential approximation. Our calculations show that, in both SmFeAsO and LaFeAsO, the Raman and IR 
phonon frequencies are well predicted by pseudo-potential density functional methods. This confirms the capability 
of ab-initio methods to describe the dynamical properties of these materials. 

We find that the phonon frequencies in the two materials are very similar at low ui, while those at high u> are much 
higher in SmFeAsO than in LaFeAsO. In particular, the E fl oxygen mode is predicted to be at 434 cm -1 in LaFeAsO, 
at 503 cm -1 in SmFeAsO. We have shown that this is related to the RE atom being closer to the oxygen plane in 
SmFeAsO than in LaFeAsO. As a general property, the same effect can be anticipated in other oxypnictides. In fact, 
results on other i?£FeAsO compounds reported in the literature show that the same mode is observed at 450 cm" 1 
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for i?i?=Ce[3lj], at 485 cm 1 for i?j?=Nd[3^|. Bearing in mind that the superconducting transition temperature T c of 
the related F- doped systems at optimum doping depends on the RE element (T c = 25 K for RE = Lap}, T c = 40 K 
for RE = Ce0,[33, T c = 51 K for RE = NdpL and T c = 55 K for RE = Sm 2]) and that the phonon frequencies are 
not strongly affected by a low level F-doping[32| , the increase of T c appears to be correlated to the increase of the E 5 
phonon frequency. Although the electron-phonon coupling constant is very low for LaFeAsOflsT]. and probably for all 
oxypicnitides, the above comparison suggests that a role of the phonons in the pairing mechanism cannot be simply 
ruled out, especially in view of the complex and strong magneto-phonon interaction [33( expected in these materials. 
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FIG. 1: Raman spectra collected between 150 and 250 cm -1 in different sample points (P1-P4) at 300 K. The 100 K spectrum 
collected at point P4 is shown for comparison. Vertical dashed lines mark the room temperature peak frequency values. In the 
inset, a typical room temperature Raman spectrum of SmFeAsO up to 600 cm -1 is shown. 




FIG. 2: (Color on line) Optical conductivity at selected temperatures in the phonon region. The 100 K peak frequency values 
are reported. The inset shows the reflectivity R(w) of SmFeAsO in the whole frequency range at the same temperatures 
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SmFeAsO SmFeAsO LaFoAsO LaFeAsO Involved atoms 

Exp. Calculated Calculated Exp. [9, 19} Simmetry Polarization Activity (RE — Sm , La) 
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TABLE I: Room temperature peak frequencies in cm - for the Raman and infrared phonons of SmFeAsO, Present room 
temperature Raman data coincide with those of Ref. [l9| . where also the Raman line at 345 cm^ 1 (marked by *) is reported. 
Peak frequencies at 100 K are given in brackets. Calculated pho non frequencies are given for SmFeAsO and LaFeAsO, the latter 
compared with literature experimental values from Refs. [9j, [l9j. Mode symmetry, polarization (X,X: in-plane, Z: out-of-plane) , 
optical activity (IR or Raman) and dominant atoms involved in each mode are also given. 



